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基于 自 适 应 控制 的 八 个 混沌 系统 的 多 级 组 合同 步 


IE, F WB. IRE 
(郑州 轻工业 学 院 电气 信息 工程 学 院 ， 郑 州 450002) 


摘 要 : 针对 传统 的 混沌 系统 结构 样式 单一 和 系统 变量 少 的 问题 ， 设 计 了 一 种 混沌 电路 系统 。 该 电路 有 较为 复杂 的 动 
力学 行为 ， 较 高 的 敏感 性 和 较 强 的 抗 干扰 性 。 随 着 电路 参数 的 变化 ， 结 合 自 适应 稳定 性 判 据 和 混沌 运动 理论 详细 分 析 
讨论 了 该 系统 的 不 同 动力 学 行为 ; 此外， 根据 电路 图 搭建 了 相应 的 电路 并 利用 示波器 观察 该 电路 的 动力 学 行为 ， 其 行 
为 与 MATLAB 的 仿真 结果 相 一 致 ， 进 一 步 证 明了 电路 的 可 行 性 和 灵活 性 。 在 此 基础 上 ， 文 章 重 点 提出 并 研究 了 基于 
本 系统 的 多 级 组 合同 步 ， 通 过 构造 不 同 的 控制 器 ， 实 现 了 八 个 系统 之 间 的 多 级 组 合同 步 。 仿 真 结果 表明 ， 该 同步 方案 
在 收 和 你 速度 和 精度 上 具有 很 好 的 效果 。 
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Multi-stage combination synchronization of six chaotic systems based on adaptive control 


Sun Junwei, Li Nan, Wang Yanfeng 
(College of Electrical & Information Engineering, Zhengzhou University of Light Industry, Zhengzhou 450002, China) 


Abstract: Aiming at the problem of single structure style and less system variable for traditional chaotic system, we proposed 
a chaotic circuit. The circuit has more complex dynamic behaviors, higher sensitivity to variable and stronger anti-interference. 
With the change of circuit parameters, we discussed the different dynamic behaviors of the system in detail based on the 
adaptive stability criterion and chaos motion theory. In addition, we constructed the corresponding circuit and observed the 
system dynamics phenomena by oscilloscope. The dynamics phenomena are in accordance with the simulation results of 
Matlab, which further proves the feasibility and flexibility of the circuit. On this basis, we put forward and studied the 
multi-stage combination synchronization based on this system. By constructing different controllers, we realized the 
multi-stage combination synchronization among the eight systems. The simulation results show that the synchronization 
scheme has a good effect on convergence speed and precision. 


Key words: chaotic circuit; period; bifurcation; transient chaos; combination synchronization; system error; controller 
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是 从 电路 设计 的 角度 来 构造 电路 ， 相 对 比较 复杂 ， 因 此 有 必要 


引言 


从 实践 方面 去 研究 和 构造 一 些 易 于 实现 的 电路 。 


随 着 非 线性 控制 的 发 展 ， 混 沌 系统 的 控制 和 


应 用 成 为 非 线 另 一 方面 ， 将 混沌 同步 应 用 到 保密 通信 成 为 近年 来 保密 通 


a 


学 的 一 个 重要 研究 方 


， 尤 其 是 以 混沌 同步 为 基础 的 保密 。 信 技 术 的 研究 热点 和 竞争 较为 激烈 的 研究 领域 。 保 密 通 信 


通信 
非 线 
同 

数 有 


伦 效 


tF 


混沌 系统 的 动力 行为 和 特征 也 不 尽 相 同 。 常 见 的 非 线 性 函 ERITEM, 如 伪 随 机 性 、 对 初始 条 件 和 参数 的 高 敏感 性 、 


[10, 11 


在 实际 应 用 中 得 到 了 广泛 的 关注 。 在 非 线性 动力 系统 中 ， 的 主要 目标 是 使 通信 系统 中 所 传输 的 有 用 信和 号 不 容易 被 窍 取 且 
性 函数 会 引起 系统 的 混沌 行为 。 如 果 系 统 的 非 线性 函数 不 “能够 在 通信 系统 的 接收 端 有 效 地 恢复 出 来 。 混 沌 信号 具有 许多 


平方 函数 、 立 方 函 数 、 分 段 线性 函数 、 单 函数 等 ， 非 线性 ，” 谱 性 和 长 期 的 不 可 预测 性 。 这 些 特 性 满足 了 一 些 通信 系统 对 通 


多 混沌 系统 的 设计 方法 都 是 基于 一 些 著名 的 混沌 系统 ， 如 洛 


Ey 
党 
E: 


有 Geneiso RA, Sprott RA”, Henon RA”, 信保 密 信号 的 特殊 要 求 , 因此 , TREE [8] 2o 7 ERR AA TR 
Arneodo-Coullet £5 £i", Chua Zi", Elwakil 系统 中 等 。 然而 ， 有 潜在 的 应 用 前 景 。 
混沌 同步 主要 是 指 响应 系统 的 状态 趋 近 于 驱动 


x 


系统 ”、Réssler 系统 ”、Chen 系统 ”等 。 这 些 系 统一 般 都 。 系统 的 状态 。 目 前 ， 常 见 的 混沌 同步 包括 完全 同步 ”、 相 位 同 


收 稿 日 期 : 2018-05-13; 修 回 日 期 : 2018-07-18 ”基金 项 目 : 国家 自然 科学 基金 资助 项 目 (61572446, 61472372, 61472372, 61603348) ; 河南 省 科技 
创新 人 才 计 划 资 助 项 目 (174200510012) ; 河南 省 高 等 学 校 重点 科研 项 目 〈162300410220) 
作者 简介 : 孙 军 伟 〈1984-) ， 男 ， 副 教授 ， 博 士 ， 主 要 研究 方向 为 非 线性 系统 的 控制 与 建 模 . 生物 信息 处 理 与 控制 ; 李楠 (190872 ， 男 ， 硕 士 研究 生 ， 


主要 研究 方向 为 非 线 性 系统 的 控制 与 建 模 、 电 力 电子 ; 王 


延 峰 (1972-), 男 ( 通 信 作 者 ) ,教授 ,博士 ,主要 研究 方向 为 生物 信息 处 理 与 控制 (810920179@qq.com) . 


201811.00165v1 


chinaXiv 


4p. REER DU XqHDIEU T. yaa JOBBIIbCU 
等 。 然 而 ， 大 多 数 的 研究 都 集中 在 单 驱 动 系统 和 单 响应 系统 的 
同步 方案 上 ， 研 究 多 驱动 系统 和 多 响应 系统 以 及 各 系统 之 间 的 
卓 合 同步 ”内 容 较 少 ， 因 此 有 必要 对 多 系统 方面 的 研究 给 
更 多 的 关注 。 

本 文 设计 了 一 种 混沌 电路 系统 ， 该 电路 有 较为 复杂 的 运动 
轨迹 ， 较 高 的 敏感 性 和 较 强 的 抗 干扰 性 。 同 时 ， 重 点 研究 了 基 
于 本 系统 的 多 阶段 组 合同 步 ， 通 过 构造 不 同 的 控制 器 ， 实 现 了 
八 个 系统 之 间 的 多 阶段 组 合同 步 。 实 验 验 证 表明 了 该 系统 有 着 


Ne 


T 


复杂 的 动力 学 特性 ， 所 提出 的 同步 方案 具有 很 好 的 准确 性 和 可 
行 性 。 

1 ”一 种 新 的 混沌 系统 及 其 基本 特性 

1.1 电路 模型 


混沌 系统 由 下 面 的 三 阶 方程 式 来 表示 。 


m 
»—b(o-x-cy)tyz 
ż=c-y 

其 中 : a,b,c 为 系统 (1) 的 参数 ，x，y，z 是 系统 (1) 
的 状态 变量 。 当 a=1, b=1, c=1 时 ， 系 统 (1) 的 Lyapunov 指 
数 如 图 1 所 示 ， 可 以 看 出 Lyapunov 指数 为 LE1>0，LE2<0 和 
LE3<0， 表 明 此 时 系统 CIO 处 于 混沌 状态 。 当 初始 条 件 : 
x(0)—1L y(0)—1L z(0)-—l1HN[, RAR (1) 的 混沌 吸引 子 
如 图 2 所 示 。 


25 


李 雅 普 诺 夫 指 数 


-0.5 1 1 1 1 1 1 1 1 1 
ü 10 20 30 40 50 60 70 80 90 100 


时 间 CO 
图 1 系统 式 (1) 的 Lyapunov 指数 
Fig.1 Lyapunov index of formula (1) 

1.20 系统 参数 的 影响 
在 这 一 部 分 运用 分 岔 图 、Lyapunoyv 指数 等 方法 详细 讨论 了 
RAR O) 的 参数 对 系统 的 动力 学 行为 的 影响 。 
1.2.1 参数 a 的 影响 

当 aE[0，5]， 参 数 b=1，c=1， 初 始 条 件 x(O) = 1L 
y(O) =1, z(O) 王 工 时 的 Lyapunov 指数 图 和 相应 的 分 倪 图 如 
图 3 和 4 所 示 。 

观察 图 3 和 4 可 以 看 出 ， 当 参数 ae [0，0.7] 时 ， 最 大 的 
Lyapunov 指数 LE1 > 0, 而且 LE1 从 正 值 渐 渐变 为 零 , 如 图 (5)(a) 
所 示 ， 当 a=0.2 时 ， 系 统 式 〈1) 产生 了 混沌 现象 。 当 参数 4 继 
续 在 区 间 [0.7,5] 增 加 时 ， 系 统 再 次 进入 混沌 状态 ， 通 过 在 区 间 
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[0.7,5] 任 意 取 值 ， 可 以 观察 到 系统 的 混沌 状态 ， 如 图 5(b) 所 示 。 
为 了 研究 初始 条 件 对 系统 式 (1) 的 影响 , 在 图 5 中 分 别 用 红色 
1 线 代 表 初 始 条 件 为 (0.8, 1, 1.4) 时 的 运动 轨迹 ， 蓝 色 曲 线 代表 
初始 条 件 为 (1,1,1) 时 的 运动 轨迹 。 


z 平 面 


(c) (d) 
图 2 系统 式 (1) 的 混沌 吸引 子 


Fig.2 Chaotic attractor of formula (1) 
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参数 a 
图 3 参数 a 的 Lyapunov 指数 


Fig.3 Lyapunov index of parameter a 
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参数 a 
图 4 参数 a 的 分 命 图 


Fig.4 Bifurcation diagram of parameter a 
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(a)a=0.2 
图 5 参数 a 变化 时 系统 式 (1) 的 混沌 吸引 子 

Fig.5 Chaotic attractor of formula (1) when parameter a changes 
1.2.2 参数 b 的 影响 

当 be[0，5]， 参 数 a=1，c=1， 初 始 条 件 x(0) — 1, 
y(O) — L z(O) 天 工时 的 Lyapunov 指数 图 和 相应 的 分 侈 图 如 
图 6 和 7 所 示 。 

混沌 系统 式 (1) 随 参数 b 的 变化 产生 了 丰富 而 复杂 的 的 动力 
学 行为 ， 如 图 6、7 所 示 。 观察 图 6 可 知 ， 当 参数 b 在 区 间 不 断 
增加 时 ， 最 大 的 Lyapunov 指数 LE1 从 正 值 变 为 零 ， 其 动力 学 
行为 由 混沌 变 成 周期 性 行为 。 而 当 参 数 be [0.1,0.2] 和 b € [2.6， 
2.9] 时 ， 此 时 最 大 的 Lyapunov 指数 LE1 < 0， 系 统 式 (1) 出 现 
了 周期 倍增 现象 ， 如 图 8(a) 所 示 。 系 统 式 〈1) 的 混沌 行为 主要 
集中 在 当 参 数 be [0,0.1] U[0.2，2.6] U[2.9，5] 这 三 个 区 间 内 。 
因此 ， 系 统 的 共存 分 岔 行为 和 共存 吸引 子 主要 发 生 在 € [0, 0.1] 
U[0.2，2.6] U[2.9，5] 这 三 个 区 间 内 ， 如 图 8(b) 所 示 。 这 里 ， 
分 别 用 红色 曲线 代表 系统 1) 在 初始 条 件 为 (1.5, 2.1,1.3) 时 系 
统 式 (1) 的 运动 轨迹 ， 蓝 色 曲 线 代表 初始 条 件 为 (1,1,1) 时 的 运 
动 轨迹 。 


(b) a=5 
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图 6 参数 b 的 Lyapunov 指数 


Fig.6 Lyapunov index of parameter b 
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参数 b 
图 7 参数 b 的 分 贫 图 


Fig.7 Bifurcation diagram of parameter b 
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(a)b-0.15 
图 8 参数 b 变化 时 系统 式 (1) 的 混沌 吸引 子 

Fig.8 Chaotic attractor of formula (1) when parameter b changes 
1.2.3 参数 c 的 影响 
M ce [0.7，5$]， 参 数 c=1，b=1， 初 始 条 件 x(0) = 1L 
y(O) — L z(0) = 1 WF} Lyapunov 指数 图 和 相应 的 分 贫 
图 如 图 9 和 10 所 示 。 
通过 分 析 图 9 和 10 可 以 看 出 系统 最 大 的 Lyapunov 指数 
LEI 从 零 变 为 正 值 ， 表 明 系 统 在 起 点 处 处 于 暂 态 阶段 。 随 着 参 
数 c 在 区 间 [0.7，1.4] 中 增加 ， 系 统 式 〈1) 进入 了 混沌 状态 
始终 保持 ， 如 图 (11) (@O 所 示 。 当 参数 ce [1.4，1.9] 时 ， 如 图 9 
所 示 ， 系 统 式 CO 的 最 大 Lyapunov 指数 LE1 等 于 零 ， 此 时 系 
统 处 于 周期 倍增 阶段 ， 如 图 (11) (b) 所 示 。 当 ce[1.9，5] 时 ， 最 
X Lyapunov 指数 LE1 为 正 值 ， 系 统 式 (1) 再 次 进入 了 混沌 状 
So o ÆR 11 中 分 别 用 红色 曲线 代表 初始 条 件 为 (0.5, 1, 2.2) 时 
的 运动 轨迹 ， 蓝 色 曲 线 代表 初始 条 件 为 (1,1,1) 时 的 运动 轨迹 。 


(b) b=1.5 


1 1 1 
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3 
参数 c 
图 9 参数 c 的 Lyapunov 指数 


Fig.9 Lyapunov index of parameter c 


图 10 参数 c 的 分 岔 图 


Fig.10 Bifurcation diagram of parameter c 
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z 平 面 


(a)c=1.2 
图 


(b) c=1.5 
参数 c 变化 时 系统 式 (1) 的 混沌 吸引 子 


Fig.11 Chaotic attractor of formula (1) when parameter c changes 
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2 ”电路 的 具体 实现 和 结论 


通过 利用 


沌 系统 式 COD 


区 


S| s 


Multisim 构造 系统 式 〈1) 的 
的 动力 学 特性 


电路 
示波器 观察 实验 现象 。 如 


Va 


实现 了 泥 


12 所 示 , 混沌 电路 由 电阻 、 电 容 、 运算 放大 器 和 乘法 器 等 构成 。 
" m " 
Ri wd Te me B 
A m R10 " 
R4 U3 | 4 | 二 
iin | J YS r lil " 
全 R7 b us 3 d ij - f | 
XY 9 | [ E. 
m 5 
P 55-5 21. z 
I Eh 35 os 
. PU. D 
x nu RS " 
图 12 Multisim 电路 图 
Fig.12 Multisim circuit diagram 
其 中 ， 运 算 放 大 器 、 电 阻 器 和 电容 器 及 其 外 围 电路 可 用 于 
实现 积分 环节 ; 运算 放大 器 ， 电 阻 器 和 外 围 电路 用 来 实现 比例 
运算 和 加 法 运算 ， 通 过 乘法 器 实现 了 各 信号 之 间 的 乘法 运算 。 


有 a=1，b=1，c=1， 因 此 系统 方 


根据 系统 aa CD 的 方程 ， 
FE (1) 可 以 构造 为 


LI R, 
x=—— 3 —— y 
1000R,C,R, 

» 1 R. R. R. 

" 5 i =i ; 5 2 
?ONWORGR CR, RO gs Qo 

, EX — x? 

1000R,,R,C, 

对 比 式 (1) 和 (2) ， 保 持 相 对 应 的 系数 相等 ， 因 此 可 以 


得 到 如 下 结论 。 


R, EN R pu R pen 
1000R, R,C, 1000R, R,C, 1000R,R,C, ” 

R; E Ro -- - AN 口 | 、 
1000R, R,C, 1000R,R,C, .: 电阻 的 阻 值 分 别 为 
R-R,-R,-100KO, R, =R, =R, =R, =R, 

R, =R, =Ro =R =R =R =10KQ ; E 7 值 是 
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根据 所 选 的 器 件 ,根据 Multisim 
。 然 后 利 


第 37 卷 第 1 


[| 


BER 
j 示 波 器 调试 , 得 到 了 系统 


名 


13 所 示 。 


[| 


很 明显 ， 示 波 器 所 得 到 的 


p Zi 
TR X 


14 5j Matlab 


2 相 一 致 ， 证 明了 混沌 系统 式 〈1) 的 正确 性 和 可 行 性 。 


vr)aod [B/A [a] [Ac] o |oc 


Channel B 
Reverse 
Save 
Ext. trigger 
Channel A Channel B. Trigger 
Scale: |1 VD Scale: |1 VDiv Ede: [r]X| [A] 
Y pos. (Div): |0 Y pos.(Div): |0 Level: [o E 
CAE Snge | Normal | Auto 
(a) x 轴 与 y 轴 的 相位 图 


(a) Phase diagram of x-axis and y-axis 


< 


Ti 图 国 | samaA 


T2 e+ 
T2-T1 
Timebase 
Range:|5 5 
X pos. (Div): 

Y/T Add B/A 


(b) Phase diagram of x-axis and z-axis 


Channel B 

E Reverse 
Save 

Ext. trigger 
Channel A Channel B Trigger 

Scale: |1 V/Div Scale: |1 V/Div Edge: [F]* [A] 
Y pos. (Div): [0 Y pos. (Div): [0 leve: [o [V 
AB] [ac] o oc [ac] o pc ||- Single | Normal | Auto [None | 
(b)x 轴 与 z 轴 的 相位 图 


< 


| Channel A Channel B 
|T: (6 E z Reverse 
Tes 
T Save 
T2Ti Ext. trigger 
| Timebase Channel A Channel B Trigger 
Range:|5 s | Scale: |1 V/Div | Scale: |1 V/Div Edge: t [4] 
X pos.(Div): 0 Y pos. (Div): |0 Y pos.(Div): |0 Level: fo y 
| |v/r|Add| B/A |Afe| |Ac| o Dc AC| 0 DC - Single | Normal | Auto |None 


(Oy 轴 与 z 轴 的 相位 图 


(c) Phase diagram of y-axis and z-axis 


图 13 示波器 调试 图 


Fig.13 Debugging diagram of oscilloscope 
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车 于 稳定 性 理论 ， 考 虑 以 下 多 个 带 控制 器 的 混沌 系统 
Ni 三 (二 机 (3) 
X2 = f,(x,) + u, (4) 
yy = g(y)*wv (5) 
Ya 7 £80) t; (6) 
az-h(z)tw (7) 
z2 =h,(z,)+w, (8) 
y, 01) *$ (9) 
W =0,(W,)+ő, (10) 
其 中 ， XX, Yo Yo Za Z Wis Y ER” , 
X = (Xi X1 X3 Xn) X = (Xs X33, X3 7X4) 
y, 2 Qi? Yu ya om Y5 = (Yas Y2» Y23 Ys) 
T T 
z= (Zi Zio Zia Zin) s Z2 = (Eis Z229 Z23 734) 
V, = Qus Vii Vis Wn), V; = (Wais VW Wan) 


分 别 是 系统 (3) — C80 的 状态 向 量 。 
fis fas 8182, 15,01, 0, 分 别 是 连续 函数 ， 
U, Ua, Vis Vas Wis Wos 所 ,多分 别 是 要 设计 的 控制 器 。 

定义 1 分 别 考虑 系统 式 (3)(4) 组 合 和 系统 式 (5)(6) 组 合 的 
误差 ; 系统 式 (5)(6) 组 合 和 系统 式 (7)(8) 组 合 的 误差 ; 系统 式 (7)(8) 
组 合 和 系统 式 (9)(10) 组 合 的 误差 ， 系 统 式 (9)(10) 组 合 和 系统 式 
(3)(4) 组 合 的 误差 ， 即 


e -(6,6€,,€,,€,). 


(11) 
@ = (65,64, 65 vena) e, 一 (epy ez ea * €, £z 

e, = (€13, ez 633 ` “63) ,es 一 (ea z4 E34 ** eu) , 
e, = (Kyo, + Kx, mici 

ez 一 (yy tb; yz —myZz; —m2221)， 

€4; 一 (M,Z; 十 mn2rZ2 —pyV; 一 P2w2z) s 

€4; 一 (pv; + Ponai — kx — kx) ， 

C Ky, Kais li; bi, mi. mo. Pis pu; 分 别 是 常 系数 ， 


i 二 1, 2, 3 )。 如 果 时 间 t 趋 于 无 穷 时 ，Jlimllel| = 0 ， 那 么 系 


统 式 (3)(4) 组 合 和 系统 式 (5)(6) 组 合 达到 同步 ， 系 统 式 (5)(6) 组 合 
和 系统 式 (7)(8) 组 合 实现 同步 ; 系统 式 (7)(8) 组 合 和 系统 式 (9)(10) 
组 合 形成 同步 ， 系 统 式 (9)(10) 组 合 和 系统 式 (3)(4) 组 合 实现 同 
步 。 


为 了 更 加 具体 地 说 明 问 题 ， 以 系统 式 (1) 为 例 并 且 考 虑 下 
且 八 个 带 控制 器 的 三 维系 统 。 
Ži = X2 + Ul 
Ži =X +X Xo + Ui (12) 


RI 2 
ži =1— xh + Us 


. 
X21 = X, t UA 
. 


X2 = —JX, X5 F XX, F Uy 


(13) 


2 
xz =l- x5 + uy, 


Yu = ya 十 VH 


yi 三 一 1 二 yao 十 yaoya3 十 mo2 (14) 


2 
Yi =1— yp 十 Vs 


Ya 7 Yn + Vo 


(15) 


Ya 5Yana F Yn F yy)23 + Vo; 


. 2 
Y3 =l- y5 + V3 


Zi —Zj 十 Wi 
Z12 = Zi E Zi F ZZi + Wiz 


2 
Z13 三 1 一 Z2 + ws 


(16) 


Z2 = Z22 + Wl 


Z2 = =Z * Za 十 Z22Z23 十 Wy, 


(17) 


: 2 
Z23 =l- zj + W3 


Va 三 Wi 十 与 | 


Wi 三 一 Vi 十 Wia VAS + G2 


Vis -1-w, 十 es 


(18) 


V3 7V35tós 


V» =W Wan Yt e» (19) 


V» -l-w; + 623 


Hr]: Us Uais Vj, Vajs Wio Wais Gi, Ġo; (1 二 1,2,3) 是 
需要 设计 的 控制 器 。 
八 个 混沌 系统 分 别 考 虑 下 面 的 多 阶段 组 合同 步 控 制 器 。 


U, z Kyu, + kb 一 Lv 一 lzivi 

U, lv, t Lv, - muw, — mw, (20) 
U, = M, Wi +M, Wz -piós -Paón 

U, = Prigi + Pag 一 Ki — Kat; 

U; = Ktn t Kw -l vi; —122 V2 

U; = liva lv; — MW — M, W3 Q1) 
U, = Mp Wo + mw, — pio, poo» 

Us = Poi + paó» — Kit; Kt 


录用 定稿 IE, F: 基于 自 适 应 控制 的 八 个 混沌 系统 的 多 级 组 合同 步 


U, = Kus 十 Ku, 一 Lavs 一 lv 
Uio —lv, 十 Lava — 3Wis — M33 W33 
U,, = maw + mw — P1313 — P23623 


U,, = p63 + pos — Kitts — K45 


根据 系统 (120 — 19) ， 多 阶段 同步 控制 器 可 以 设计 为 


(22) 


以 下 形式 ”23-25 
U, 2 -k,x — Kx» +hiy + ly yy 7 06, + Be 
U, =- Y2 Lys t mz; + mz, - Per 0,6, Q3) 
U, = -Mi Man tp tp - Be — fa 
U, = Pu 7 Pal» +K, +k, + Qe * Bes Tey 
U; kx, — KX — KX tk Xo — Ky Xo 
= kx» - liy + liz Yiz ths yi 13 — 1 yo1 
+ 1 Ya +ly Ys Y23 7 0565 + Y;65 
U, 2ljy, -loyo -loYoYi + lys 7l» —122 22 23 
Z M21 + M212 M2212213 ~ M3222 + M327 
Hm zz 一 De (24) 
U, = M321 — M2 7 M2212213 M221 ~ My27 
710525 25, 7 PV FPW HPY nA 
7 Pa3V 1 + P23V FPW 25V 53 — /3€3 
Us = PW n HPV» PioVioV is 7 PV» * PV 
+ pal; 5; — Ki Kj + Ki 7 K5xj 
KX — KoX5X — 06564) — V261 
U,- -ka+kuax -katkat tla -ly 
+l az 70565 + fies, 十 13643 
Uio = -la+lay7 1 一 La+la3y n Tm, 一 maz 
+M; 一 maz 一 ea = Yen 
Q5) 


2 2 
U 2-m;-4 mz -m, +M, » +p1s — 
2 2 
Pi Pt P3 7 PY n + 06564 = Y365 
2 2 
U, = —Pis + PY ,-PatPsV ntks a 


2 2 
Kx, + ky = KX » 70563 — 0564 — 565 


RP, 05.8, 7 CSL 2,3) 是 正常 系数 。 

定理 1 如 果 系 统 误差 式 (11) 的 组 合 控制 器 
Us FUE T IN Us UU Us UU ,UU 
空 制 器 满足 条 件 式 (23) ~ (25) ， 那 么 系统 式 (12) ~ (19 
将 达到 多 阶段 组 合同 步 。 


证 明 ”构造 如 下 的 Lyapunov 函数 : 


vod gu 


4 
Va)- 22. (e2  e2 +e2) (26) 
i=4 


V7) 对 时 间 t 求 导 可 得 


. D . . . . 


. 
V(t) 2 €, eit € ente; ete €21 65, 0m 65, ent 
4 ; . š ; n 


€, 631+ € € €, 633 十 £4, 04 Heyy €32+ €, 633 (27) 


考虑 系统 式 〈12) ~ 〈19) 可 以 得 到 以 下 结论 : 
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ei = Ky xit Kk, Xa od yi 一 Di Ya 


= ky + ky =l Yn 一 Di 
+ ku + kun, 一 iv =l, v 


$ E . : 
en zl, yy tl yj 7 my zu- ms za 


=l yp tl y» 7 mz, 7 M722 
liv t b; = mw, 7 mw; 


è 
€3i =M, Zu +M; Z217 pi Vn- Pa Ya 
= M + mz puli Pal» 


MW +M, Wy 一 DG 一 De 


ea 2 py Ynt py V 7 Ki Xi k, Xni 
= puo + pai Kx — Kx» 
Pui + P2621 一 KU Kit 


en 2 Kj xit ky x2=l yp 一 y» 
= -kX + kX Kx; 7 kun 
tk, tksxyxs + ly 一 ay 
=la Yag + by 7155 ly)» 
+ ki + kd 一 ly 一 7 


en =1 Yntly y2 一 mi zo- m, zn 
=la Yn +h + hy, yy 7 liy, 
+y Ya *15y5y,*miz, 7 mz, 
一 mxzpza * IZ, 710525 = Illy zs 


lv + lv = mW, — mw» 


€» — my zit m, zn- py Yr Py V» 
= mz, + M222 + M2212213 7 M3221 
+My +M HPW 7 Po 
SPW * PV ~ Pr 一 P22W22723 


+m Wy +My Wy 7 Phi -Pasa 


en 2 py Ynt Pa V7 Kj X2- Ky, Xn 
= poV 7 Poi» “PW os * P»V ~ PV» 
7 pU Xs Kx 7 Kx, Kta 十 Ki 


nn 7 Ky XXy3+p12612 + Pan 7 Killa 7 Kt; 


@13 = Kj X13+ Kj, x23- li Yla y 
2 2 2 
= kg -Kat ks -kux » =la lay — 
2 
l3 thy stk Kus lav; 7 lav 


en zl yl y47 m zo- m; zs 


2 2 2 2 
=la l3) 1 +ly 一 3 2 一 ma+mazi 一 mo 


2 
+ mz 2 十 1373 Ll = miaWs — M, W23 


€33 = Dy zit m5 257 pi V i57 Pos V a3 
2 2 2 
= M Z M37 12 + M3 = M3 22 = Pis PY n27 Pa 


2 
十 Da Hm Wi, + mW 7 pis 7 Dass 


es = pi Vit Pos V 557 Ki 357 ky X23 
2 2 2 
= pi = PW i t Pa Pa 5 7 Ks t Kx 5— 


2 
kj + kx; +M,3W,3 + msw, - kun — kt; 


Q8) 


Q9) 


G0) 


通过 式 (23) ~ (250 对 式 (28) ~ (300 简化 ， 可 以 得 到 
如 下 的 结论 : 


201811.00165v1 


chinaXiv 


2 2 2 2 2 
=V 3) 7 05€ 45 7 05€ 44 — Pye 34 7 Y3€ 05 


2 2 
一 ]3e 3705€ 43 


<0 
于 ,Bi, Vi (i 二 1,2,3) EEX RZ, VO 是 非 负 


数 ,所 以 yG) 是 非 正 值 ,根据 Lyapunov 稳定 判 据 , 系统 误差 (11) 


随时 间 趋 于 无 穷 时 ，lim||e| = 0 ， 证 明 完 毕 。 


4 “仿真 结果 


仿真 实验 证 明了 所 提出 的 理论 的 正确 性 。 假 设 
a =l, J, =2, 7 =3, œ, =2, B, 21, y, =4, 
a, =3, B =4, 7 =1, k; =h; =m, = P; = 
k, =L; =m; = Pa =1,(i=1,2,3)。 系 统 的 初 
台 状态 任意 给 出 如 下 所 示 : 
(x1 C0), x,,(0), x3 (0) =, 2.4) ， 
Qx3, (0), x5, (0), x23 (0))" = (2,5, D”, 
GO). yO, y (0) = (2, 4,3)", 
(0), y0), y (0) =2,5, D", 
(zı (0), z2 (0), z,(00))" = (2, 4, 3)”, 
(z4,C0), z,,(0), z,,00))" = (6, 2,3)”, 
C0), v5 (0). v44(0))! = 4, 5, D7, 
(Wa C). Wa (0), Wo3 (0) = (4, 2,537, 
根据 式 (9) ， 可 以 计算 出 系统 误差 的 初始 值 分 别 是 
(ei1(0), e,,(0). ez, (0), 6,00)" = 1, — 4,3, 2)”, 
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(ea(0), e, (0), e, (0, 64,0)" = (-2, 3, 一 1,0)7， 
n0; 1€21 
. f (ei3(0), e54(0), e€44(0), e44(0))" -— (2. —2, —1,2) ,系统 
ena =- e-e a LL 
: uiid G1) 的 误差 收敛 图 如 图 14 所 示 。 
63 = = fe T7183 
es = Qez + Be — 6 
€i = —056, F yes 
en = —D,€, 
i 2 pb 22 (32) 
€32 = —y,€, 
en — —05€4 — Nen 
À 40 2 4 ; à 10 12 16 48 20 
en 2 -Ae + Be t Ys 时 间 (s) 
€ 7 — fe, — ye, 
. 3513 ^ /3€3 (33) 
633 = 0364 — 3€ 
€43 == 0463, 05€, — 36s 
将 式 (31)~(33) 分 别 代 入 式 (27) 可 以 得 到 
V(t) 2 e, eut €; e» 64 €5 €), 21 €, 00 65, 5 
Ez 6E31 十 €4, €32 t C33 €33 £4, C41+ £4, € 42 + €, €43 f ， ， 
2 2 2 2 2 ?8 1 2 4 5 6 8 9 10 
=- Qe 17e 3 7/4 Qe bp 一 De 22 (34) PHI) 


图 14 


误差 收敛 图 


Fig.14 Error convergence diagram 


5 ”结束 语 


本 文 设计 了 一 种 新 的 混沌 电路 ， 通 过 自 适应 稳定 判 据 理论 
详细 讨论 了 系统 参数 对 于 系统 动力 学 行为 的 影响 ， 显 示 出 系统 
具有 较 好 的 混沌 特性 。 其 次 ， 重 点 提出 并 证 明了 一 种 新 的 同步 
方案 一 多 级 组 合同 步 , 通过 构造 不 同 的 控制 器 , 结合 误差 分 析 ， 
实现 了 八 个 系统 之 间 的 多 级 组 合同 步 。 实 验 结果 表明 该 方案 具 
有 一 定 的 可 行 性 。 
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